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In his recent work on a tilt instability for advanced LIGO interferometers, P. Savov discovered 
numerically a unique duality relation between the eigenspectra of paraxial optical cavities with non- 
spherical mirrors: he found a one-to-one mapping between eigenstates and eigenvalues of cavities 
deviating from flat mirrors by h{f) and cavities deviating from concentric mirrors by —h{r), where 
h need not be a small perturbation. In this paper, we analytically prove and generalize this re- 
markable result. We then illustrate its application to interferometric gravitational-wave detectors; 
in particular, we employ it to confirm the numerical results of Savov and Vyatchanin for the impact 
of optical-pressure torques on LIGO's Fabry-Perot arm cavities (i.e. the tilt instability), when the 
mirrors are designed to support beams with rather flat intensity profiles over the mirror surfaces. 
This unique mapping might be very useful in future studies of alternative optical designs for LIGO 
interferometers, when an important feature is the intensity distribution on the cavity optics. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

Laser Interferometer Gravitational-wave Observatory 
(LIGO) and other long baseline detectors, are formed 
by high-Finesse Fabry-Perot arms in order to increase 
the circulating optical power and to enhance sensitivi- 
ties by suppressing shot noise. LIGO intcrfcromcrom- 
eters, as well as the baseline design for advanced-LIGO 
detectors (whose funding is scheduled to begin in FY 
2008) 3], all use spherical mirrors. Non-spherical mir- 
rors have been alternatively studied and are considered 
for use in gravitational-wave interferometers, for the sup- 
pression of thermal noise they offer. 

In particular, as shown by O'Shaughnessy, Thorne and 
Agresti 0, IE IE S 13 > the thermal fluctuations of mirror 
surfaces are better averaged over by a flat beam, such 
a meso-like Q mode. The corresponding optical design 
has shown a strong tilt instability ^3 ^.nd Thorne has 
proposed a different version of the mesa beam, that is sup- 
ported by nearly concentric and opportunely shaped mir- 
rors; this new version provides the same intensity profile 
at the cavity mirrors (and thus the same thermal noise), 
but imply a weaker tilt instability (even weaker than cav- 
ities with nearly concentric spherical mirrors analyzed by 
Sigg and Sidles 0, 0|) — as calculated by Savov and 
Vyatchanin ■ A general method to design a family 
of optical cavities has been devised by Bondarescu and 
Thorne jl3j||, from nearly flat resonators to nearly con- 
centric ones. 

Mesa beams are constructed by coherently overlapping 
Gaussian beams, with either(i) translated parallel axes, 
or (ii) axes in different directions but sharing a common 
midpoint Ts?]. Mirror shapes which support such beams 
as fundamental modes are derived from the phase fronts 
at the mirror locations, with case (i) corresponding to 



Mexican-hat mirrors, and case (ii) corresponding to the 
nearly-concentric version. Using the resulting optics pro- 
file, higher-order optical modes and eigenfrequencics of 
the designed cavities must be calculated by solving an 
eigenvalue problem, which has been done for nearly-flat 
cavities by O'Shaughnessy and Thorne 0,0,|EI3|j and for 
nearly-concentric cavities by Savov and Vyatchanin [lOj | . 
During his numerical work, Savov discovered that the 
deviation of nearly-concentric Mexican-hat mirrors from 
concentric surfaces is exactly the opposite of the devi- 
ation of nearly-flat Mexican-hat mirrors from flat sur- 
faces; he also found that the corresponding higher modes 
of these cavities all have the same intensity profiles, and 
that there is a one-to-one mapping between their eigen- 
values. He went on and conjectured a general duality 
relation between axisymmetric cavities with two identi- 
cal mirrors facing each other: cavities with mirrors de- 
viating by — ft,(|r I) from concentric surfaces (nearly con- 
centric mirrors) will support modes with the same in- 
tensity profiles and related eigenvalues as cavities with 
mirrors deviating by /i(|r |) from flat surfaces (nearly flat 
mirrors). It should be noted that the deviation /i(r) is 
not required to be inflnitesimal, it can change the mirror 
shape arbitrarily as long as the paraxial approximation 
is still satisfied^. While such a duality relation is well- 
known between cavities with spherical mirrors, i.e., those 
with ft,(r) ar^ (for example see llJ, HE EE 113 ) ; to our 
best knowledge no such relations had been established 
between generic cavities. 

In this paper, we prove this remarkable correspondence 



^ Here and henceforth in the paper a 2-D vector f has been used to 
indicate each point on planes orthogonal to the the cavity axis. 
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analytically, for a even broader category of cavities: those 
whose mirror shapes remain invariant under the parity 
operation, identified as spatial reflection in the two di- 
mensional r-space (which is also equivalent to a f80° 
rotation around the cavity axis). Eigenmodes of such 
cavities can be put into eigenstates of parity, and we 
show that all corresponding eigenmodes of dual cavities 
have the same intensity profiles at the mirrors, with their 
eigenvalues satisfying 

^fc^(_l)P. + lg-2.feL(^fc)*^ (1.1) 

where (— 1)^''= is the parity of the fcth eigenmode; sub- 
scripts c and f denote nearly concentric and nearly flat 
mirrors, respectively. 

We will give two alternative proofs of this duality re- 
lation. The first one relies on the geometrical properties 
of the propagator from mirror to mirror. In this descrip- 
tion the eigenfunctions are field amplitudes at mirror sur- 
faces, and we see right away that the corresponding eigen- 
states have the same intensity profiles there. The second 
proof is based on the "center-to-center" propagator. The 
center-of-the-cavity fields are the eigenstates and the cor- 
respondence relation is manifested by a two-dimensional 
Fourier transform, that univocally relates the dual cavi- 
ties. 

This paper is organized as follows. In Sec. ^ we report 
the first proof; in Sec. Ill Al the Cartesian coordinates are 
used and some general features of the eigenproblem are 
described; in Sec. Ill Bl the cylindrical coordinates are 
used, and the case of axisymmetric resonators is stud- 
ied. Section IIIII contains the second proof and the 2- 
D Fourier transform relation between the center-of-the- 
cavity eigenmodes of dual cavities. Section llVl specializes 
to the case of Mexican-hat cavities. When the nearly- 
flat and the nearly-concentric mirrors are implemented 
in the system, the corresponding mesa beams are con- 
nected by Fourier transform, as we report in Sec. IIV Al 
In Sec. IIV ill plots and analytical forms are provided, for 
the amplitude distributions at the center of the cavity 
and at the mirror surfaces; in Sec. IIV CI we address the 
tilt instability of the nearly concentric Mexican-hat res- 
onator and show how easily it can be analyzed, applying 
the duality relation to the results obtained for the nearly 
flat Mexican- hat cavities 1^ - We comment and re- 
view the implications of the general duality in Sec. 

II. ANALYTICAL PROOF FOR 
MIRROR-TO-MIRROR PROPAGATION 

A. In the Cartesian coordinate system 

In this section we focus on field distributions on mirror 
surfaces, and restrict ourselves to cavities with two iden- 
tical mirrors facing each other. The extension to the non 
non symmetric cavity is presented in Appendix A. We 
adopt the Fresnel-Kirchoff diffraction formula to prop- 
agate fields from mirror surface to mirror surface (see 



e.g. 23). In this formalism, the field amplitude wi(r') 
on the surface of mirror 1 propagates into 

V2(r) =y"dV'/C(f,f')wi(r') (2.1) 

on mirror 2, via the propagator 

IC{r,r') = ^il + cos0)e~'''P k='^, (2.2) 

Airp A 

from r ' (on mirror 1) to r (on mirror 2), where p denotes 
the (3-D) spatial distance between these two points and 

6 stands for the angle between the cavity axis and the 
reference straight line, as is illustrated in Fig. ^ We 
know that the Fresnel-Kirchoff integral eigenequation 

jv{r) = J d^f' JC{r,r')v{r') (2.3) 

univocally determines the eigenmodes v and eigenvalues 

7 of the cavity. 

h(r') 




FIG. 1; Symmetric Nearly Flat Mirrors. 
Applying the paraxial approximation 

e^O, + '^p' - h{f) - h{r') , (2.4) 

and we can use 

IC'^(r,r') = ^e-»fcie'feM?)e-i^l-^--^'l'e*'=''(-"') . (2.5) 
2TrL 

in the integral eigenequation. 

Here the mirror surfaces deviate by h{f) from a flat 
reference, and the subscript f is used to reflect this con- 
vention. From here on, we will also refer to /C[' as the 
nearly flat propagator. We now consider two slightly de- 
formed concentric mirrors (see Fig.|2J) so that the mirrors 
height with respect to the flat reference surface is 

h{f')=f^/L +b{f'), (2.6) 

where the height b{r) is the deviation from the concentric 
spherical surface (note that concentric spherical mirrors 
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have their radii of curvature equal to i/2, and thus sur- 
face height jV). Inserting Eq. (|2.6|1 into Eq. (|2.4() . we 

obtain the propagator for a nearly- concentric cavity, 



K\{rA) 



— ikL 



ik 

gikb{r)^+H\r+f'\^^ikb(r') 



(2.7) 



We use the term nearly concentric propagator for 
/C^(r, r'). Although we use the terms nearly-flat and 
nearly- concentric, h and b are not required to be small; 
in fact, they can represent any deviation from perfectly 
flat and concentric spherical mirrors. 
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FIG. 2: Symmetric Nearly Concentric Mirrors. 

Now let us consider mirrors that are then invariant 
under parity, i.e., those in which we also have 

h(f) = /i(-r) , h[r) = 6(-r) . (2.8) 

so that /Cf_c are both invariant under a spatial reflection 



{r, r '} <-> {— r, —r'^ 
and therefore, we have 

VK. = K.V , 

where we have defined 

Vv{f) = v{—r) . 



(2.9) 



(2.10) 



(2.11) 



for two dimensional reflection. Equation H2.10|l implies 
that all eigenmodes can be put into forms with definite 
parity. We derive the following relation between nearly 
flat and nearly concentric propagators, as constructed: 



[lC^i^r,r')]A-e 
that is equivalent to: 



2ikL V — h 



K-\ry), 



(2.12) 



(2.13) 
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Nearly Concentric 
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Eigenstate 


Vt 


Vi 


Parity 


AY 


AY 


Half-trip eigenvalue 


7f 




Round-trip eigenvalue 







TABLE I: Correspondence of propagation kernels, eigen- 
states, parities, and eigenvalues between dual configurations. 

Suppose we have an eigenstate Wf of /C^, i.e., an eigen- 
state of a cavity with mirror deviating by from flat 
surface, and we compute its eigenvalue 7t and know the 
parity eigenvalue (—1)^: 



7f Df , 

Afvi- 



(2.14) 
(2.15) 



By applying Eqs. 
correspondance 

AAt 



^J^-^J^, we derive the 

g-2ifcL(_-^)p+l^*^*^ (2.16) 



which identifies Vc = Vf as the corresponding eigen- 
state of IC~'^, that is eigenstate of the corresponding 
resonator we denote the dual. The eigenvalue is 7c = 
e~^*'^^(— l)^'+^7f . We also induce that the parity is still 
(— 1)P. The reverse is straightforward and the result is 
an established one-to-one correspondence between dual 
cavities. We summarize this mapping in Table ^ It is 
obvious to note that that the corresponding eigenstates, 
V[ and Vf , have the same intensity profiles on the mirror 
surfaces; for infinite mirrors, we know V{{f) is real- valued 
(see Appendix 01, so it is an eigenstate of the dual con- 
figuration itself. 

For cavities with identical mirrors facing each other, 
the full, round-trip propagator is just the square of the 
half-trip one. From Eqs. H2.12(l and H2.1()|l . we have 



[/C: 



„4ikL 



Ac 



(2.17) 



which means that the same duality correspondence ex- 
ists between eigenstates of the full propagator, with their 
eigenvalues related by 



Vc 



(2.18) 



Note that when h{r) = r^/(2L) the two dual cavities 
are identical to each other. Using the relation that links 
the eigenvalues of two dual resonators, we can determine 
the spectrum 



7c 



-ikL-\'in7T /2 



where n G Af. The resulting separation between the 
eigenvalues is the Gouy phase 



gi^G arccos(l — L/i?) 



R = L 

computed for confocal resonators [TsllTGLIr^ . 
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B. Specializing to cylindrical mirrors 



The mapping of the eigenvahies are given by 



In most LIGO appHcations, cavity mirrors still have 
cylindrical shapes: h{r) = h{\f\). This allows us to 
decouple radial and azimuthal degrees of freedom, and 
simplify the eigenvalue problem. We shall follow roughly 
the notation of 

We adopt the cylindrical coordinate system: 



r = r(cos ip, sin (p) 



(2.19) 



Since JC is now invariant under rotation along the z-axis, 
all eigenmodes can be put into eigenstates of rotation: 



v{r, if) = R{r)e , m = integer . 



(2.20) 



Inserting this into the eigenequation ^2.1^ and performing 
analytically the angular integration we obtain the radial 
eigenequation 



K\^){r,r')Rnm{r'ydr' , (2.21) 



where for each angular mode number m we have indexed 
the radial eigenstates by n, and 



f(m)(''j "f"') — 



L 



' Jrr 



krr' 



^jfc[-L+ft(r)+/«(r')-i%^] 
(2.22) 

is a symmetric radial propagator, in the nearly-flat de- 
scription.^ Since K^^^^{r,r') is symmetric, we obtain or- 
thogonality relations between radial eigenfunctions: 



Rmm{r)Rn2'm{r)rdr = S^r 



(2.23) 



Comparing Eqs. H2.24|l and (|2.22|l . we obtain: 



(-1) 



m+l 



At 



(m) 



^ 2^kLJ^-h 

c(m) 



(2.24) 



(2.25) 



This is a radial version of Eq. (|2.13(l : here we know ex- 
plicitly that all TO-eigenstates have parity (—1)™. 

Following a similar reasoning as done in the previous 
section, Jor each angular mode number m, we can estab- 
lish a one-to-one correspondence between radial eigen- 
states of a nearly-flat configuration to those of the dual 
configuration: 



[R-n 



[R-n 



(2.26) 



[in 



(-1) 



m+l -2ikL 



[in 



(2.27) 



Similarly, the round-trip eigenstates have the same cor- 
respondence, their eigenvalues related by 



[Vn 



^ L'/nmJf 



(2.28) 



III. ANALYTICAL PROOF BASED ON 
CENTER-TO-CENTER PROPAGATION 

A. Propagators for vacuum and mirror surfaces 

In this section, we focus on complex amplitudes of the 
optical field on planes perpendicular to the optical axis 
(the z axis). An optical mode propagating along one di- 
rection of the optical axis can be specified completely by 
the distribution of the field on the z = const plane. For 
example, we denote the optical field on the plane z = zi 
by v{f,zi), where f is the 2-D coordinate of the point on 
this plane. The effect of any linear paraxial optical sys- 
tem (including open space, thin lenses and mirrors) with 
input plane zi and output plane Z2 can be characterized 
by its transfer operator, U, which takes the form of an 
integration kernel: 



w(f,Z2)= / dr'U{r,Z2;r',zi)v{r',zi). 



(3.1) 



In particular, the operator that describes the paraxial 
propagation down a length L in vacuum is 



Using Eq. H2.6|l again, for a configuration with b{r) 
correction from concentric spherical mirrors, we obtain 
the radial kernel of the nearly- concentric description: 



GL{r,r') = i 



2ttL 



e''-''^ exp 



—ik 



(f 



2L 



(3.2) 



For a mode propagating in the ±z direction with field 
(complex) amplitude distribution v{r',zi) at z = zi, the 
amplitude distribution on a surface described by height 
z(r ) — zi ^ ^(^) is given by 



v[r, zif')] 



(3.3) 



Here we emphasize that the spatial point of interest is 
located outside the z = zi plane, and that the 2-D vector 
r describes the projection of that point onto the z = zi 
plane. 

From Eq. 1)3.3(1 . one deduces that the operator for re- 
flection off a perfect infinite mirror with shape /i(r) is"^ 



T^[h{f)]{f,r') = -S{r - f") 



-fl\ 2ikh{r) 



(3.4) 



It is easy to verify that both Ql and TZ\^h(r)] are unitary 
operators. 



2 Here we have used Jr,(z) = — '— /^^ e'^ '^e*"^cZip, where 

27ri" Jo 



3 



The minus sign in Eq. 13. 4i is used because we use a convention 



Jn(z) is the rath order Bessel function of the first kind. in which a phase shift by tt is gained upon reflection. 
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B. Analytical proof based on center-to-center 
propagation 

In this section we present an alternative proof moti- 
vated from the construction of the flat-topped beams 0, 
Eo| : (i) the nearly flat configuration has its fundamental 
mode generated by spatial translation of minimal Gaus- 
sian beams, while (ii) the nearly concentric configuration 
is generated by rotation (of propagation direction at the 
center of cavity) of minimal Gaussian beams, or a trans- 
lation in the momentum fc -space. This had led us to 
speculate that the two sets of eigenstates correspond to 
each other via Fourier transform (similar to the relation 
between position and momentum space in quantum me- 
chanics) . 

We will use the operator Gl/2 [see Eq. (|3.2() ] which 
propagates the field forward by half the cavity length. 
For simplicity we denote it by fj: 



g{r,r') = z^e-*'=^/2e- 



■ik 



(3.5) 



Using g and Tlh{r) [defined in Eq. (|3.4|l . with h{r) the 
mirror surface height], we can re-express the eigenvalue 
problem as: 



^[h{r)]U = G'R-[h(r)\Qu = 7U , 



(3.6) 



with >C[/j(^)] the center-to-center propagator when the 
mirror deviates from flat surfaces by /i(r ), in which the 
optical mode propagates from the cavity center to the 
mirror, gets reflected, and propagates back to the center. 
In fact, C is related to the surface-to-surface propagator 
/C by a unitary transformation. 



(3.7) 



This means the two proofs are mathematically equiva- 
lent. Similar to )C, the operator £ also commutes with 
parity or [Cf. Eq. fTTUIl ] 



(3.8) 



With the propagator on hand, we proceed with our in- 
tuition that the modes must be related by Fourier trans- 
forms. In order to do so, we first define the 2-D Fourier- 
transform operator as 



•^(r ,r ) = ^ 

which satisfies 

It is easy to show that, 
[g*T-^] (r,f') 



(3.9) 



(3.10) 



7r2i2 

= [^e''^-^7^[^.2/(2L)]e] (r,r'). 



(^2-"g¥[r "-(r-r")]' 



(3.11) 



[The integral on the second line can be done by inserting 
a factor e"'''' ^ into the integrand, and then letting e — > 
0+.] Similarly, [or by taking the transpose of Eq. (|3.11() ]. 
we have 



^Q*ie^'^^Q7i[p2/(^2L)] ■ 
Using Eqs. (|3.11() and H3.12|l . we have 



(3.12) 



2ikL 



-e 



[rV(2i)]'^[-'iA]'^[?V(2i)]i^ 

T-'C^n^^T. (3.13) 



Here and are mirror heights related by the duality 
relation. 



hA{?) + hB{?) = r''lL, 
and we have used the fact that 



(3.14) 



7?.[r V(2L)]^[-/iA]^[r V(2i)] = ^[rVi-''A] = "^['ib] ■ 

(3.15) 

According to Eq. H3.13|l . given any eigenstate ua of 
^[Ha] with eigenvalue 7a and a definite parity of p, we 
have 

(-1) 7a«a = >'^[hA]^A 



2ikL 



^Cyu^^{Tu\) = (-1) 
In other words, the mapping 

UA ub — ^u\ 



T-'Cyn^}{Tu*A), (3.16) 
f+i--2^fe^7X(J^u^) . (3.17) 



(3.18) 



transforms each eigenstate of 'Cf/iA] i^ito its dual one of 
£[/i3]; the corresponding eigenvalue relation is 



7B = (-l)^+^e-2*^-^7l. 



(3.19) 



For similar reasons, given any eigenstate ub of U\ji^-^ 
(with definite parity) , Tu'^ must also be an eigenstate of 
U[h^]- Moreover, since 



(3.20) 



the state Tu^ is in fact the inverse image of ub [under 
the mapping H3.18|l ]. This means we have established a 
one-to-one correspondence between eigenstates of C[hA\ 
and those of >C[/ib]- 

Now let us look at intensity profiles on the end mirrors 
surface. For the eigenstate wa, the field amplitude at the 
constant- 2; plane of the end mirror is Qup^. For its image 
eigenstate ub = ^u\, we have 

guB = g{TuX) = [g*T-'uA]* 

= [ie'''''R[r2/(^2L)]5uA]* (3.21) 



6 



which does have the same intensity profile [see Eq. (|3.4|) ]. 

For the round-trip propagator C? ^ using Eqs. (|3.13|) 
and H3.8|) . we have 



Here a is the waist size, which we leave general (rather 
than setting a = cTmin) for the moment. The duality 
image of Vi is 



^2 



AikL 



1 pi '77 



(3.22) v,{r) = {Tvl\{r) 



so we have the same duality correspondence (I3.18f) be- 
tween eigenstates of the full propagator, with the map- 
ping between eigenvalues given by 



(3.23) 



IV. APPLICATION OF THE DUALITY 
RELATION USING MESA BEAMS AND 
MEXICAN-HAT CAVITIES 



d foe ^ T 



roEV 



roev 



with 



L 
2k 



2'Ka 



r2. 

mm 5 



(4.2) 



(4.3) 



The mesa beams were constructed to have flat-topped 
intensity profiles at the cavity mirrors with rapid fall- 
off near mirror rims, in order to achieve lower thermal 
noises 0, IE IE ■ There are two versions of mesa beams 
with the same intensity profile, the nearly flat and the 
nearly concentric. Cavities that support them (Mexican- 
Hat cavities) are related by the duality relation, as real- 
ized by Savov , during his study of radiation-pressure- 
induced tilt instabilities. In this section, we shall explic- 
itly construct these two fundamental modes, study their 
relations at the center of the cavity, and at the cavity 
mirrors. We will also discuss analytical features of the 
two modes that have not been obtained before. We will 
also give an example of how the calculation of the tilt 
instability can be dramatically simplified for nearly con- 
centric Mexican-hat cavities, using the duality relation, 
based on results already obtained for the nearly fiat con- 
figuration. 



When going from Eq. (|4.1|l to Eq. H4.2() , the Fourier trans- 
form has been completed by two steps. First, the spa- 
tial translation by r o is replaced by the phase factor of 
e L ^ which represents a tilt of the propagation axis 
by an angle of 2r o/^- Second, the cr-Gaussians turn into 
(T:„-Gaussians. [This correspondence between Gaussians 
in fact reflects the duality between pairs of spherical cav- 
ities.] As a consequence, Vc represents the superposition 
of Gaussians with symmetry axes going through the cav- 
ity center, but with tilt angles distributed uniformly in 
a disk with radius 2p/L — exactly the construction of 
a nearly-concentric mesa beam. In particular, Eq. I|4.3|) 
tells us that minimal Gaussian would have turned into 
itself in this process. Hence we have shown explicitly 
the correspondence between the nearly-flat and nearly- 
concentric mesa beams (the fundamental modes of the 
corresponding cavities) . 



A. Construction of Mesa beams in Cartesian 
coordinate system 

Nearly-flat Mesa beams are constructed by coher- 
ently superimposing minimal Gaussians, namely Gaus- 
sian modes with the sm allest p ossible spot size at the 
cavity mirrors, (Tmin = \/ L/ {2k), whose axes are parallel 
to the cavity axis and lie within a cylinder centered at 
the cavity axis. At the middle of the cavity, the axes 
intercept with the constant-z plane in a disk V, with ra- 
dius p. It is evident that such a construction will give 
a rather fiat intensity profile in the central region of the 
end mirror with radius ~ p; beyond this radius, the in- 
tensity profile falls off as a Gaussian with decay length 
(Tmin, which is conceivably the fastest possible IEIEQ- 

The complex amplitude of the nearly-fiat mesa beam 
(fundamental mode of the corresponding cavity) at the 
center of the cavity is of the form 



B. Profiles of mesa beams and mirror shapes 

In order to study mesa beams in more details, we adopt 
the cylindrical polar coordinate system (r, 0); the cylin- 
drical symmetry of these beams will make the complex 
amplitude only depend on r. Equations 14.1|l and (|4.2|1 . 
written in the polar coordinate system, become 



rodro 



"0 "'0 
2tt 



r — 2rorcos((/)-0o)-|-r'Q 
TTT^ 



"0 "'0 



r 4-22ror cos((/) — (/)o) 

d<po e 



(4.4) 



(4.5) 



Vf{f) 



roev 



1 



2tt(t 



(4.1) 



Here wq = \/L/k = v^cTmin and L is the total length of 
the cavity. Carrying out the angular integrations analyt- 
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ically, we get 



p/tuo 



2a;oe-(="'+^«)/o(2xa:o)da;o, (4.6) 



Ji{2xp/wo) 



(4.7) 



On the other hand, if r is much larger than wq [since 
p ^ wq, this region overlaps with the previous one], we 
can apply the asymptotic expansion of Iq 



(4.9) 



where x = r/wo, and Iq is the modified Bessel function 
of the first kind. Examples of normalized power distri- 
butions of nearly fiat and nearly concentric mesa beams 
are plotted in the upper panels of Fig. O In these plots, 
we take p = Awq , which corresponds to the configuration 
proposed for Advanced LIGO (for reasons that will be 
explained in Sec. lIVC!|l . 
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on Eq. (|4.6|l . and obtain 
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-(p/wQ-x) 



,(4.10) 



where we have defined x = r/wQ. From Eq. H4.10|l . we 
note that when wq ^ r <C p, we recover the result of 
1; when r gets close to p, the amplitude will drop, 



..waist 



similar to the tail of an error function. Qualitatively, 
we could write Wf-Mcsa(p) ^ P- In the ultimate limit of 
p/wq +00, we have 



Vf 



^{r) = 1, p/wq -> +00 . 



(4.11) 



The concentric configuration, on the other hand, has 
a completely different field distribution. According to 
the analytic expression H4.7|l . the amplitude must be dis- 
tributed within a radius of x wg/p 1, or r '--^ ^o/Pi 
which is much smaller than the waist size of the minimal 
Gaussian. In this case, we could also qualitatively write 
Wc-Mcsa(p) ~ Wq/p- In the limit of p ^ 00, we use 



Ji{ax) 



S{x) , a — + +00 



(4.12) 



and have 



FIG. 3: Comparison between nearly flat (left panels) and 
nearly concentric (right panels) Mesa beams. Upper panels; 
normalized intensity profiles at the center of the cavity. Mid- 
dle panels: normalized intensity proflles at mirror surfaces 
Lower panels: phase fronts at the position of the mirrors. 

Let us analyze these amplitude distributions in more 
details, in the case of p 3> wq, i.e., when we translate the 
minimal Gaussians by a distance substantially greater 
than their waist size. For the nearly-flat configuration, 
we can easily see from Eq. (|4.1ll that, when (p — r) / wq S> 
1, the field distribution can be approximated as 



^(r) = S{x) , p/wo-* +c 



The fact that 



Wf-Mosa(p) • fi'c-Mosa(p) ~ 



(4.13) 



(4.14) 



clearly refiects the Fourier-transform relation between 
two Mesa beams with the same p. 

Now, let us turn to field distributions at the cavity mir- 
rors. Applying the propagator between parallel planes in 
the polar coordinate systems (eq. (j^J), 



Wf (r < p) « / cPfo 



1 



27rCT 



(4.8) 



G(r',0';r,0) 

-ikL/2 -ik\r^+r'^~2rr' cos(<p' - 4>)] / L 

■kL 



(4.15) 
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we obtain the fields 

Jo Jo 

_ r l+i 1 r r'^-2ror' cos(0'-0o)+rg 1 
e ^0 \, (4.16) 

Jo Jo 

r r r'^+2iror' cos(0' ~0o)-zrp 1 

e ^0 J , (4.17) 

at distance L/2 from the waist. Comparing Eqs. (|4.1t)|) 
and H4.17|l . we have 

[v^'^ir)] * = e'''=^^'/^v™'^(r ) . (4.18) 

It is then obvious that the two beams have the same 
intensity profiles on the cavity mirrors: 

\vr''ir)\^\vT''{r)\. (4.19) 

(An approximate formula for the end-mirror intensity 
profile was given in the Appendix of 5].) We plot these 
intensity profiles at the mirror surfaces in the middle pan- 
els of Fig.El 

Let us now determine mirror shapes by imposing that 
the optical phase is constant (which we take as for 
simplicity) on each mirror surface. We have 

y^nd(^^ykh,ir} ^ \v^''<^{r)\^ (4.20) 
^end^^ykh^r) ^ (4 21) 

Taking the complex conjugate of Eq. 14.20|l . and combine 
with Eq. (|4.21|) . using Eqs. (|4.18|) and H4.19|l . we have 

htif)^^-K{f), (4.22) 

which is the duality relation between mirror surfaces. In 
the lower panels of Fig. 21 we plot the shapes of mirror 
surfaces, again, we assume p = Awq. 



corr. from flat corr. from cone. 



-0 .25 
-0 . 75 
-1 . 25 




FIG. 4: Flat mesa beam wave front (left panel) with respect 
to a flat surface and concentric mesa beam wave front (right 
panel) with respect to a concentric surface, as analytically 
computed. 



C. Applications of Mesa beams to Advanced LIGO 

In order to achieve lower thermal noise in the test 
masses, the intensity profiles at the mirrors must be as 
fiat as possible. In the case of infinite mirrors, the choice 
is to use cavities with flat or concentric spherical mirrors, 
whose eigenmodes have uniform (absolutely fiat) profile 
distribution. However, the mirrors must have finite sizes 
(e.g., as limited by the size of the beam tube), and the 
intensity profiles must be confined to a very large extent 
within the rims of the mirrors, in order to decrease the 
diffraction loss upon each reflection. In Advanced LIGO, 
a power loss below lOppm is required 0. For this reason, 
we are forced to deviate from flat or concentric configu- 
rations — to such an extent that the diffraction loss is 
within the requirement. When only spherical mirrors are 
used, if on the one hand we decrease the radius of cur- 
vature from -|-oo (fiat), and on the other hand increase 
the radius of curvature from L/2 (concentric), the dual 
configurations, with 

l/(2i?i)-f l/(2i?2) = (4.23) 

will have the same intensity profiles at the end mirrors, 
thus the same diffraction loss and thermal noise. For ex- 
ample, i?i = 54 km and R2 = 2.077 km both give exactly 
the loss specification, while i?i is the current baseline de- 
sign. However, spherical cavities are not optimal in terms 
of their thermal noise: (the two types of) mesa beams, 
whose intensity profiles are ffatter given the same loss 
specification, turn out to provide much lower thermal 
noises |^ For these beams, the larger the param- 
eter p, the lower the thermal noises, but the higher the 
diffraction loss. The loss specification of Advanced LIGO 
corresponds to p = 4wo which is the case we study in 
Fig. El 

While having the same diffraction losses and thermal 
noises, dual configurations do differ significantly in a 
very important aspect — their eigenspectra are different. 
Thus, any problem using modal analysis of optical cavi- 
ties will reveal these differences and probably the duality 
relation if nearly fiat and nearly concentric configurations 
are compared. 

One such problem is the radiation-pressure-induced tilt 
instability: as the mirrors tilt, the beam inside the cavity 
walks away from the center of the mirrors, producing a 
torque, which in some cases can drive more tilt in the 
same direction, and become destabilizing (see Fig. [SJ. 
As shown by Sigg while for all cavities there is 

always one tilt mode in which the radiation-pressure- 
induced torque is destabilizing, the instability is much 
weaker in nearly concentric configurations than in nearly 
fiat ones. The reason is that while in the two cases the in- 
tensity profiles are identical, the optical axis of the beam 
walks away by a much smaller distance in the concentric 
case, given the same amount of tilt in the unstable mode 
(see Fig. jSJ . According to Sigg's calculation for spherical 
mirrors, the tilt instability for a nearly fiat configura- 
tion with Advanced-LIGO power (^ 1 MW circulating in 
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FIG. 5; Comparison of tilt instability of nearly flat and nearly 
concentric symmetric optical cavities. For more details see 

Ref. mm. 



the cavity) can be too strong to handle for the angular 
control system. For this reason, we would prefer nearly 
concentric cavities. 

For general, non-spherical cavities, a perturbative pre- 
scription for calculating the tilt instability has been for- 
mulated by Savov and Vyatchanin |lOl |. in which the tilt 
instability growth time is expressed in terms of eigenval- 
ues and intensity profiles of the cavities' spatial eigen- 
modes (Eqs. 2.13, 2.14, and 4.8 of O). Savov and Vy- 
atchanin applied their prescription to both nearly flat 
and nearly concentric Mexican-Hat cavities; in partic- 
ular, they had to solve the eigenvalue problem for the 
nearly concentric cavities in order to obtain the eigen- 
values and intensity profiles. Savov discovered the du- 
ality relation in this process. Had the duality relation 
been known, one could have taken the eigenvalues and in- 
tensity profiles of nearly flat Mexican-hat cavities, avail- 
able from previous works, applied the duality transfor- 
mation, and obtained the tilt instability for nearly con- 
centric Mexican-Hat cavities without having to solve the 
eigenvalue problem again (see Section VI of T^). 

Finally, let us make a qualitative comment on the nu- 
merical magnitudes of tilt instabilities in the various con- 
figurations considered. Numerically, according to Savov 
and Vyatchanin |lOl |. we have 

nearly flat MH (p = 4wo) 

i 

nearly flat spherical {R — 54 km) 
I ^ 

nearly concentric spherical (R — 2.077 km) 

i 

nearly concentric MH (p = Awq) 

with configurations less and less unstable from top to bot- 
tom. Interestingly, this sequence of decreasing instability 



is consistent qualitatively with the corresponding mirror 
shapes: with the same amount of diffraction loss, the flat 
MH does appear more flat than the nearly flat spherical 
mirrors, while the nearly concentric Mexican-Hat mirror 
does appear closer to concentric than the nearly concen- 
tric spherical mirror. 



V. CONCLUSION 

In this paper, we provided two different analytic proofs 
for Savov's duality relation between symmetric cavities 
with mirror height h(f) measured with respect to a 
flat surface and those with mirror height —h{r) mea- 
sured with respect to a concentric spherical surface (valid 
within the paraxial approximation): the corresponding 
eigenmodes have the same intensity profile at the mirrors, 
their amplitude distribution at the center of the cavity 
is related via Fourier transform, while their eigenvalues 
are related by complex conjugation (see Table QJ. These 
two proofs are based on the mirror-to-mirror propagator, 
and the center-to-center propagator, respectively. 

We illustrated this duality relation with the two types 
of Mesa beams proposed for Advanced LIGO. In partic- 
ular, we showed explicitly that these beams are related 
to each other by a Fourier transform at the center of the 
cavities, and that they have the same intensity profiles at 
the end of the cavities. We also related the mirror shapes 
of the Mexican-Hat cavities that support these two modes 
by the duality relation. In addition, the duality relation 
could have allowed us to avoid solving the eigenequations 
once more for the nearly concentric Mexican-Hat cavities, 
and used instead results already available for nearly flat 
Mexican-Hat cavities. 

The duality relation can also be applied to the more 
general Mesa beams by Bondarescu and Thorne, which 
interpolate between nearly flat and nearly concentric cav- 
ities: in these beams, minimal Gaussians are both trans- 
lated and have their propagation axes rotated, to differ- 
ent extents jl3] . 
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APPENDIX A: DUALITY RELATION FOR 
NON-IDENTICAL MIRRORS 

In this section we will study the duahty relation when 
the mirrors shapes are not identical, but each still sym- 
metric under a 180° rotation around the cavity axis. 
Since now the field distributions of eigenstates over the 
two mirror surfaces are not the same, we have to study 
the eigenvalue problem associated with the round-trip 
propagator, instead of the individual surface-to-surface 
ones. Nevertheless, we can still use the propagators (|2.5f) 
and 1)2. 7|l to build a system of integral equations relating 
field distributions wi(ri) and V2{r2) over the two mirror 
surfaces. [All through this section, we use the subscripts 
1 and 2 to refer to quantities associated with mirrors 1 
and 2, respectively] If the mirrors deviate from parallel 
planes by /ii,2(?^), we have: 

7i^'i(^i)= / <i^r2 ICi2{ri,r2)v2{r2) , (Al) 

JS2 

72^2(^3)= / <i^ri IC2i{r2,ri)vi{n) , (A2) 

JSi 

where 71^2 are the "eigenvalues" and 

K:i2(ri,r2) = j^|ri-r.p+»fc/».(r.)^ (^3) 

2t:L 

2TTL 

are the propagators from mirror 2 to mirror 1, and from 
mirror 1 to mirror 2, respectively. The equations (|A1|) 
and HA2|I give the field at each mirror in terms of the 
reflected fleld at the other but they can be combined to 
form the round-trip equation which states that the field 
at each mirror must reproduce itself after one round-trip. 
In the following, we will add a subscript f or c to make a 
distinction between quantities related to the nearly-fiat 
or nearly-concentric case. 

Vfv.fin)^ J^^d\-[JC^,fHn,n)vif(f[), (A5) 

Vf V2f{f2) = J^^ dVl JC'^f^ (f2, r-) V2f(f'2), (A6) 

where the common eigenvalue »7/ is given by 71/72/ and 
the round-trip propagators 

/C^i''^(fi,f'i) = / d'^'r2 K.12 fin, 1^2)1^21 f{r 2, n) 

¥i\p{?2.r'2) = (1 - 2) . K\p{r^,T\) (A7) 

In the nearly-concentric configuration, using kernels of 
the form (|2.7|l for the propagation from one mirror to 
the other and combining them as done for the nearly-fiat 
configuration, we obtain the following nearly-concentric 
round-trip equation for the field distribution over the 



mirror 1 (similar formula for the mirror 2 with the sub- 
stitution 1 2). 

VcViM)^ f dVi/Cj^''^(fi,f;)i;ie(f;) (A8) 
ICt'Hnyi) - -/^ dV2 e-'^''{^y ■ (A9) 

. gir|f'l-(-r2p-|-^|f='2+r'iP+ifc6i(fi)+ifc6i(r="i)+2ifc62(?2) 

where 61^2 are the mirrors deviations from concentric sur- 
faces. Using the assumed symmetry properties of the 
mirrors, the propagators for the nearly-flat and nearly- 
concentric cavity fulflUs this relation (the same is true for 
the mirror 2 with the substitution 1 <-> 2) 

= e-*'''^[IC'^}''^n,f[)]* (AlO) 

Equation (|A10|) . together with Eqs. (|A7|) and HASp . pro- 
vides us with a more general duality relation, for cavities 
with non-identical mirrors: as long as the corresponding 
mirrors of two cavities A and B satisfy 

hc^Air) = - Ksir) , a = 1,2, (All) 

the eigenstates and eigenvalues of the two cavities will be 
related by: 

v^a = vIb, VA^e-^''''^7^*B, a =1,2. (A12) 

APPENDIX B: EIGENSTATES AND 
EIGENVALUES FOR CAVITIES WITH INFINITE 
MIRRORS 

When the mirrors are infinite, it is straightforward to 
check that two fundamental properties, 

y"dV'/C(f,r")/C*(r',f") = S{r~r"), (Bl) 

IC{f,f') = IC{f',r), (B2) 

are satisfied by both propagators ICf and K,^; they can 
be re-written into 

/C/C^ =1, /C = /C^ , (B3) 

where I is identity operator, iC^ the conjugate of /C, and 
/C^ its Hcrmitian conjugate. In simple terms, K. is unitary 
and symmetric. It is well known that for unitary oper- 
ators, all eigenvalues have modulus 1, and that eigen- 
vectors with different eigenvalues are orthogonal to each 
other. 

Now suppose we have an eigenvector v, with eigenvalue 
7i 77* — 1- By complex conjugating the eigenequation 
ICv — jv, we obtain 

IC*v* = -i*v* = -f-^v* ; (B4) 
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using Eqs. HB3|) . we have JC* = K,^ = JC ^, and hence 

JC-^v* = 7*w* ^ 71;* = JCv* . (B5) 
This means v* and w are both eigenvectors with eigen- 



vahie 7. We can then replace v and v* by two real eigen- 
vectors of the eigenvalue problem, v + v* and {v — v*)/i. 
This corresponds to the physical fact that the optical 
phase of eigenstates must be constant on mirror surfaces. 
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